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We show that starting with either the non-extensive Tsalhs entropy in Wang's formahsm or 
the extensive Renyi entropy, it is possible to construct the equihbrium statistical mechanics with 
non-Gibbs canonical distribution functions. The statistical mechanics with Tsallis entropy does 
not satisfy the zeroth law of thermodynamics at dynamical and statistical independence request, 
whereas the extensive Renyi statistics fulfills all requirements of equilibrium thermodynamics. The 
transformation formulas between Tsallis statistics in Wang representation and Renyi statistics are 
. . . presented. The one-particle distribution function in Renyi statistics for classical ideal gas and finite 

' particle number has a power-law tail for large momenta. 
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I. INTRODUCTION 



The Gibbs distribution function has several successful applications in various domains of physics. It is a 
consequence of statistical mechanics which satisfies all postulates of equilibrium thermodynamics 0,13 • Recently, 
however, a number of experimental data have appeared, where the asymptotic distribution differs from the 
I> [ Gibbs distribution Q • These deviations point out that either the equilibrium assumption fails for high energies 
or there is a need to construct the equilibrium statistical mechanics with non-Gibbs distribution function. 
I The Tsallis statistics explores this second alternative on the base of a non-extensive entropy j^]. There has 
been an increasing interest in studying the non-extensive Tsallis statistics. Unfortunately, the generalized 
Tsallis statistical mechanics have problems with the zeroth law of thermodynamics. Several authors 0,0 have 
, attempted to get around this difficulty, but remain unpersuasive. Further attempts to solve this problem were 
' undertaken by introducing an extensive representation for the non-extensive Tsallis entropy 0j 13 • As it will be 
shown in this work, this idea leads to a transformation of Tsallis statistics into the Renyi statistics, but does 
^ ^ not solve the problem of second law of thermodynamics for Tsallis statistics: the conventional correspondence 
I ' between temperature and heat is less obvious. In this paper we point out that the statistical mechanics with 
Oh! extensive Renyi entropy fulfill all requirements of an equilibrium thermodynamics and still have its canonical 
distribution function in a form with power-law asymptotics. It might be, on the other hand, less preferable 
regarding the stability property (H theorem based on a generalized entropy) . 

This paper is organized as follows. In sections II and III we discuss the microcanonical and canonical 
ensembles in the incomplete non-extensive statistics and in the Renyi statistics. The transformation rules from 
Wang formalism to Renyi statistics are derived in section IV. In section V the properties of thermodynamical 
Ci I averages and the zeroth law of thermodynamics are discussed. Finally, in section VI we apply these results for 
the classical ideal gas of massive particles. 

II. INCOMPLETE NON-EXTENSIVE STATISTICS 

The Wang's formalism of generalized statistical mechanics uses Tsallis' alternative definition for the equilib- 
rium entropy Q 



S ■ 

and utilizes a new norm equation |^ 0| 



kZ!0. (I) 



Here is a probability of ith microscopic state of a system and q ^ TZ (entropic index) defines a particular 
statistics. In the limit g — > 1 eq. |^ approaches the Gibbs entropy, 5*1 = — fc^p^lnpi. The q - expectation 

i 

value of an observable O in this framework is defined as follows 

= Ep?0- (3) 
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From the beginning the useful function 



should be introduced so that 



7y = ^p, = (4) 



S = k^. (5) 

1 — q 

For the derivation of the distribution functions in microcanonical and canonical ensembles we will use Jaynes' 
principle 



A. Microcanonical ensemble {E, V, N) 

In order to find the distribution function f — pi, one maximizes the Lagrange function 



$ = — - a 
k 



(6) 



For the probability pi we obtain the following expression 



p1 = 



(7) 



The parameter a has been eliminated using equation ^ and In the limit g — > 1 the distribution function 
^ has the well known formp^ = e^^^/^. Substituting eq. Q into (0) and taking into account the conservation 
rules for microcanonical ensemble results in the equality 



i-(i-,)f 



= X! ^V^,V^N,,N5E,,E = W. 



(8) 



Based on this we get the equipartition probability from eq. ((TJ as a function of the thermodynamical ensemble 
variables, the energy E, the volume V and the particle number N: 
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The entropy is given from and ® by 



S = k 



(H/i/<?)9-i _ 1 



(9) 



(10) 



Note that this entropy 1)1011 does not satisfy the famous Boltzmann principle [l^ . 

Differentiating eq. l(2Jl and using eqs. @ and (|10|l we obtain an expression for the heat, comprising the 
microcanonical limit of the second law of thermodynamics 



SQ = TdS = 0, 



(11) 



where d is the differential operator of three independent variables {E, V, N) and T is the temperature of the 
system. The first law of thermodynamics is satisfied since the heat transfer in microcanonical ensemble is equal 
to zero 



SQ = dE + pdV - fidN = 0, 



(12) 



where p and fi are the pressure and chemical potential of the system consequently. Comparing (|11|) and (|12|l 
we obtain the fundamental equation (first law) of thermodynamics 

TdS = dE+ pdV - fidN. (13) 

This equation is especially useful since the independent variables are also the environmental variables. 
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B. Canonical ensemble {T,V,N) 



The Lagrange function for the canonical ensemble is given by 



$ = a 

k 



(14) 



After maximizing the function 11411 and using eq. (^3) to eliminate the parameter a we arrive at the following 
expression for the distribution function 



[l + il-q)qf5iA-E,)] — , 



(15) 



with A = E — S/qk(3. Let us now fix the parameter /3. Differentiating the function A and eq. with respect 
to P, and using the distribution function l|15|) we obtain 



d(3 ~ ' dp- 



So the parameter /3 can be related to the temperature 

1 _ dS dS/d(3 



kp, 



T dE dE/dp 

Thus expressing P with the physical temperature T the distribution function (jl5|l takes the form 

A~E,' 



Pi 



l + (l-9) 



kTq- 



where A is determined from the normalization condition eq. Q 



5vi,v5Ni. 



N 



1 + (1 - 9) 



A-E, 



kTq- 



= 1. 



(16) 



(17) 



(18) 



(19) 



Thus it is a function of the canonical thermodynamical variables, A — A{T, V, N). Then canonical averages are 
calculated according to eq. lO in the following manner 



1 + (1 - 9) 



A-E, 
kTq-^ 



(20) 



(For another derivation of the distribution function l|18|) see the Appendix A.) 

We would like to verify now the distribution function l|18() from the point of view of equilibrium thermody- 
namics. Applying the differential operator of ensemble variables (T, V, N) on the eqs. iQ, 121 and on the energy 
E from eq. (|SJl with eigenvalues Ei, and using the distribution function l{TS|l and eqs. Q, ©, we obtain again 
the fundamental equation of thermodynamics, 



TdS ^ dE +pdV - ndN, 
where [dEi/ dT)y ]S! = and the pressure and chemical potential take the form 

'dE, 



dV 



P = Y^pI 

i 

^' = HpI (H) 



T,N 



T,V 



The free energy of the system is defined as 

F = E -TS = A + kTq-\l - ?]). 
where rj (cf . eq^J is a function of the canonical variables (T, V, N) : 

A-E, 



V ■ 



N 



l + {l-q) 



kTq- 



(21) 

(22) 
(23) 

(24) 

(25) 
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The entropy is calculated from eq. ©. Using the definition (|24(l . from eq. (|21|l we get 

dF = -SdT - pdV + iidN, 



(26) 



and some further useful relations for the entropy, pressure and chemical potential if the free energy (|24|l is 
known 



dF 

of 



-S, 



dF 
dV 



-P, 



dF 
dN 



(27) 



The energy can be calculated from the formula E = —T^{d/dT){F/T)v,N- 
Note that by introducing the new function \1/ — ^'(A) so that 

_ 1 



F = -kT- 



the Legendre transformation rule between entropy and free energy applies: 

S 1 



T— 



dT \ q-l 



(28) 



(29) 



Substituting the function A = E — q ^TS into eq. 118|) we obtain the canonical distribution function in a 
representation similar to a recent formalism by Tsallis 



Pi =V' 



1 + (1 - g) 



E-E, 



(30) 



In this case the probability l|30|l depends on two unknown variables E and rj. Thus, in order to normalize the 
distribution function, one has to solve two equations. For instance the norm equation Q, X^jP? ~ 1' 



i 

and the definition of the function = "^^Pi 



V 



SVi,vSNi,N 



E-E,, 



kTq ^rj 



E - E, 



kTq ^1] 



(31) 



(32) 



Thus E and are functions of the canonical thermodynamical variables, E = E{T,V,N) and t] — 7]{T,V, N). 
Then the entropy is calculated from formula ©, the free energy F from eq. H24() . General canonical averages 
take the form 



l + il-q)- 



E-E, 



' kTq ^f] 

or can be calculated from eqs. 127|l by differentiating the free energy F(T, V, N). 



(33) 



III. RENYI STATISTICS 

The Renyi statistics can be derived from the Renyi-entropy 

1-q 

where pi is the probability to find the system in ith microstate. The normalization is achieved by 



and the averages are taken as follows 



(34) 



(35) 



(36) 
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Here Oi is an eigenvalue of the operator belonging to the observable O. The function 

i 

and the entropy are simply related: 



5 = fc 



In 77 



(37) 



(38) 



A. Microcanonical ensemble {E, V, N) 

In order to find the microcanonical distribution function one maximizes the entropy (|34ll with an additional 
norm condition. The Lagrange function becomes 



(39) 



As a results, after using eq. (|37|l to eliminate the parameter a, the microcanonical distribution function pi takes 
the simple Gibbs form 



Pi = e 



Then from the norm equation (|35|l we obtain 

e~ = ^^Svi,vSNi,NSE,,E = W. 
i 

The distribution function H40() and the entropy are given by the familiar expressions 



and 



S = k\nW. 

The Renyi statistics in the microcanonical ensemble resembles the Boltzmann-Gibbs statistics. 



(40) 
(41) 

(42) 
(43) 



B. Canonical ensemble (T,V,N) 
The functional to be maximized in the canonical ensemble is given by 



<i> = a 

k 



After the usual procedure the canonical distribution function takes the power-law form 

E - E, 



Pi = 



kTq 



(44) 



(45) 



In this case the probability H45|l depends on two unknown variables E and 77, and, in order to normalize the 
distribution function, we solve two equations. For instance the normalization condition H35|) . J^iPi — 1^ 



r] 9-1 = 5vi,vS- 

i 

and the definition of the function rj = X^i-P? ; 



Ni.N 



l + {q-l 



E - E,, 



Ni,N 



1 + iq-l 



kTq 



E - E,, 



kTq 



(46) 



(47) 
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Thus E and 77 are functions of the canonical thermodynamical variables, E = E{T, V, N) and 77 = 7](T, V, N). 
Canonical averages (13611 take the form 

O = T]~ ^ 5vi,vSNi,N Oi 
i 

Eventually the entropy is calculated from the formula . 

Applying the differential operator with respect to the ensemble variables (T, V, N) on eq. 134(1 and using 
the distribution function H45() . and eqs. ((35|I - H37|I . we get the fundamental equation of thermodynamics 1(21(1 . 
TdS = dE + pdV — i-idN, and the condition that {dEi/dT)v,N = 0. The pressure and chemical potential are 




1 + (g - 1) 



E-E, 
kTq 



(48) 



The free energy is given by 

F = E -TS ^ E - kT^^^ (51) 
1 - 9 

and dF — —SdT — pdV + jidN . Averages can be calculated from eqs. ((27(1 and the energy from the relation 
E = —T'^{d / dT){F /T)y iq . Note that the Legendre transformations ((29(1 are valid for Renyi statistics with any 
function * = 'i'{E,r]) defined by formula (PS|) . 



IV. TRANSFORMATION RULES 



It is useful to express the demands of equivalence in the microcanonical and canonical ensembles for the 
distribution functions of Wang eqs. iQ, l(3U(l and Renyi eqs. ((42(1 . (1451) and for the averages eq. (0), ((3t)() 

[Pi r ^ Pi \ (52) 

by using the equations for the Tsallis index q and temperature T 



1 



Here the indices W and R refer consequently to Wang and Renyi statistics. Note that in the microcanonical 
ensemble {E, V, N) we obtain the following relations 

T^R) = (q(w'))-iT(^)W '^<"'' , (56) 

p(^) = pf'^), (57) 

= A^(^) (58) 



On account of eqs. ((52(1 . ((54(1 the relation between Tsallis (eq. ^| and Renyi fea. 134(1 entropies becomes 

In I ^ l«V" / — I ■ 

qi^) - 1 



l + (g(^)-l)- 



(59) 



Some authors [3, B HJI interpret this equation as the extensive representation of Tsallis entropy, but eqs. ((55(1 
and ((59() are transformation formulas from the Wang formalism of Tsallis statistics to Renyi statistics. Sub- 
stituting eqs. ((55(1 and ((59(1 into (1111) results in an invariance of the second law of thermodynamics during this 
transformation 

rWrf^W (60) 
The fundamental equation (first law) of thermodynamics eq. ((13() is invariant, too. 
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V. PROPERTIES OF THERMODYNAMICAL AVERAGES 



In this section we investigate the additivity of thcrmodynamical averages in microcanonical and canonical 
ensembles. For this reason we divide a system C = ^ + i? to two subsystems A and B and make demands of 
dynamical and statistical independence, i.e. we require that the energy of microstates are additive, 

E,,{C)^E,{A)+E,{B) (61) 

and the probability factories 

P,,{C)=pM)Pj{B). (62) 
Then the Tsallis entropy (Q) is a non-extensive variable, following rules of pseudo-additivity 0,^1 

S{C) = S{A) + S{B) - (1 - g) , (63) 

but the Renyi entropy ()34|l is an extensive variable 

S{C) = S{A) + S{B). (64) 

The average values of energy, particle number and volume (cf. 13 I36|l are extensive variables 

E{C) ^ E{A) + EiB), (65) 
V{C) ^ V{A) + V{B), (66) 
N{C) = N{A)+N{B), (67) 

under the conditions lt)l|) - (|62|l . if they follow restrictions on the number of particles and on the volume in 
the ijth microstate of the system C: Nij{C) = iV^(A) + Nj{C) and Vij{C) = Vi{A) + Vj{C), respectively. 
Substituting now eqs. (|63|l . (|55|l - Ht)7() into the thcrmodynamical relations following from eq. (|13|l in case of 
the microcanonical ensemble {E, V, N) and into the thcrmodynamical relations following from eq. H21(l for the 
canonical ensemble (T, V, N), as well as using the function r] defined in I^J, we arrive at the following relations 
between temperatures in the Wang statistics: 

TcVc = Tava = Tbtib (68) 
while the pressure and the chemical potential are equal in equilibrium: 



PC = PA ^ PB, (69) 
Mc = MA = Mb- (70) 

Here rjc = TjArjB due to statistical independence 1()2|) . The relation between temperatures in the Renyi statistics 
due to the extensive property of the entropy 164() is the conventional one: 

Tc^Ta^Tb. (71) 

The pressure and the chemical potential satisfy eq. (|69|l and eq. (|70|l . they are intensive variables. 

Thus it is proved that the non-extensive entropy of the system destroys the intensive property of the physical 
temperature T, conjugate to S. Other variables, like pressure p and chemical potential /i, conjugate to TV and V , 
remain intensive. The zeroth law of thermodynamics is violated this way (see eq. (I68II ') therefore the statistical 
mechanics with non-extensive entropy JQ) does not satisfy all principles of equilibrium thermodynamics with 
the probability pi being the equilibrium distribution function. However, Renyi statistics completely satisfies all 
these principles. 



VI. THERMODYNAMICS OF THE CLASSICAL IDEAL GAS 



Finally we would like to provide an example for the comparison of the two approaches we have discussed so 
far. Let us consider a classical non-relativistic ideal gas of N identical particles in the canonical ensemble of 
Renyi's and Wang's incomplete statistics respectively. For exact evaluation of various (/-expectation values the 
integral representation of the Euler's Gamma function will be applied, following Ref. (15f| . 



A. Wang's incomplete statistics 



The functions E and 77 are obtained by solving eqs. (|^ . for N independent particles. We get 



(72) 



and 



ri = A-\CZ,{T)y-^^, (73) 

where A = l + {1 — q)^N, Z\ is the partition function of the ideal Boltzmann gas in the canonical ensemble 0, 
[l7| . Zx = (l/7V!)(5y/A|,)^, and At ImkTYl'^ is the thermal wave length for a particle For the 

variable A we have the condition A > 0. The variable C in the case g > 1 becomes 



c 



(g(g-l))f^r(-^) 



g > 1, 



which is valid under the condition —3N/2 + q/ [q — 1) > 0. In the case g < 1 we have 

r(T^) 



c 



(q(i-g))i^r(^ + |7V)' 



g < 1, 



(74) 



(75) 



where r(z) is Euler's Gamma function. Then the entropy is calculated from eq. (5). The free energy H24|l takes 
the form 



F : 



and the pressure is obtained from eq. H23|l : 



-kT - — - + -NkTq-^ 77, 
1 — g 2 



p — — kl q n — . 



(76) 



(77) 



The A^-particle distribution function of the classical ideal gas in the canonical ensemble takes the following 
form 



/(pi,...,p^.) = -^^^ry-. 



1 ^2 



(78) 



N -n 

2m' 

which is normalized to unity 

J d^Pl...d^PN fiPi,...,PN)^l- (79) 
The one-particle distribution function is defined as an integral over the momenta of all other particles: 

/(pi) = J d^p2 . . . d^pN f{pu ■ ■ ■ ,Pn). (80) 
Integrating formula ()78|l in Wang's incomplete statistics we obtain 



f{p) = D 



\2'KmkTq ^rjA 
where the coefficient D for g > 1 is equal to 



3/2 



l-(l-g) 



2mkTq ^rjA 



A 5 
1-q 2 



D = 



(g-i)^r(^-i(A^-i)) 
r(^ - 1^) 



g> 1. 



(81) 



(82) 



(83) 



We have the restriction that —3N/2 + g/(g — 1) > 0. In the opposite case, for g < 1, we have 

p (i-g)^r(^ + fA^) 

r(^ + |(A^-i)) ' ' 

The distribution (|81|l in the — > 00 limit gives back the classical formula with an effective temperature, 
Tcs = Trj/q even at finite g — 1: 

1 \ 3/2 
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B. Renyi statistics 



The funetions E and 77 are found solving eqs. (^BJ, (|T7|) : 



E = ^TVfcT (85) 



and 

, 1-9 



ry=(^-(i+T^)CZi(T)) "\ (86) 
where A = 1 + (g — l)g^^|A^ and it should be A > 0. For the variable C in the case g > 1 we have 



C= , "^'^ , q>l, (87) 



and for one in the case g < 1 we have 



C= , ^7,, / , 9<1- (88) 

(^)^^r(^)' 

eq. ijSHIl is vahd if -3N/2 + 1/(1 - g) > 0. Then the entropy is calculated from eq. 

S'/fc = lnZi+lnC-(l + -^)lnA (89) 

1 — g 

The free energy H51II takes the form 

F = -fcTln^i + -NkT - kTlnC + (1 + Y—^)kT\nA. (90) 
and then the pressure is calculated from the second equation of (|23|l 

p^^kT^'-^. (91) 
The A^-particle distribution function of classical ideal gas takes the form 



(92) 



Integrating formula (|8U|I using eq. (|92|l for one-particle distribution function in Renyi statistics we obtain 



fip) - D 



1 



3/2 



2'KmkTqA 

where the coefficient D for q > 1 is equal 

(g-i)ir(^ + fiv ) 
r(i^ + i(^-i)) 

and in opposite case for 5 < 1 is 



2mkTqA 



qA 5 

q-l 2 



(93) 



^ = — - — 5 — ^ , a < 1, (95) 

which is valid in the hmit -37V/2 -I- q/{l - q) > 0. 

In this respect it should be noted that in the paper of A. Lavango [l^ the one-particle distribution function 
does not correspond to the iV-particle distribution function of non-extensive thermostatistics. 

All averages of Renyi statistics, eqs. (P ^ - (|?^ . can be obtained from the averages eqs. l(7^ ~ (p^ by using the 
transformation formulas of section IV. 

The one-particle distributions for Wang and Renyi statistics both approach the Gibbs distribution in the 

— > 00 limit even for q ^ I. This means that non- interacting or short-range interacting systems, like the ideal 
gas, can show non-Gibbs statistics only with finite particle number. This complies with the general physics 
experience. 

Fig.l shows the one-particle distribution function for different values of the entropic index q in the above 
described classical ideal gas of Renyi statistics. The changes in behavior of the one-particle distribution function 
as the number of particles grow can be seen in Fig. 2. 
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VII. CONCLUSIONS 

In conclusion we summarize the main idea of this paper. Two power-law tailed representations for non-Gibbs 
canonical distribution functions pi = p{Ej) were obtained, the Wang and Renyi statistics. The equilibrium dis- 
tribution function pi both for Wang and Renyi statistics satisfies the fundamental equation of thermodynamics. 
The transformation between them preserves the first and second laws of thermodynamics. The same is true 
for the Legendre transformation structure and for the thermodynamical relations for intensive quantities, like 
pressure and chemical potential in the microcanonical and in the canonical ensemble. 

As a general rule the averages associated with Renyi statistics coincide with those for Wang statistics, but 
with a renormalized temperature and entropy. This makes the physical interpretation of these quantities more 
uncertain than it was in the Gibbs case. The physical temperature in Wang statistics reflects the transformation 
formula between Wang and Renyi statistics. The non-extensive entropy of Wang formalism changes, however, 
an important property of temperature, the satisfaction of the zeroth law of thermodynamics at dynamical and 
statistical independence request. The extensive Renyi statistics, on the contrary, passes even this test fulfilling 
all the requirements of equilibrium thermodynamics. The canonical distribution function also in this case is 
power-law tailed. 

We have studied the Renyi statistics for classical ideal gas of massive nucleons with a finite number of particles. 
The equation of state (EOS) satisfies the ideal gas law {pV — NkT) familiar from Gibbs statistics in this case, 
but the distribution function differs from the Maxwell-Boltzmann one: it has a power-law tail. This difference 
may be seen experimentally only at large momenta in the tails of distributions, unless the Tsallis index is very 
near to the lower limit of gcr = 1 — 2/(3iV) for a classical ideal gas system consisting of N particles. 

Acknovifledgments: This work has been supported by the Hungarian National Research Fund OTKA 
(T034269) and by the Hungarian-Russian Inter-Academy Collaboration, as well as partially supported by the 
Moldavian-U.S. Bilateral Grants Program, CRDF Projects MP2-3025 and MP2-3045. We acknowledge valuable 
remarks and fruitful discussions with V.D. Toneev and P. Levai. One of authors, A.S.P, is grateful to the 
Department of Theoretical Physics of KFKI for the warm hospitality during his stay. 

APPENDIX A: DISTRIBUTION FUNCTION OF CANONICAL ENSEMBLE IN INCOMPLETE 

STATISTICS 

Other definition for Lagrange function of canonical ensemble is 



After maximization of (|A1|I we have 



p'i = [q{l-{l-q)PEi)] — , (A2) 
To normalize this distribution function it is necessary to introduce two new parameters instead of (3 

p'l^^[l-{\-q)(3'E,]^^, iki) 
where the parameter /?' ^ (3. Then the parameter Z is founded from the norm eq. |2l 

Z = ^ [l-(l-g)/?'i?,]T^. (A4) 

i 

Let us to find the parameter /?'. Substituting ljA3|l into eq. ^ we have 

9-1 

Z~ - 1 ,9-1 

S^-k + kl3'Z— E. (A5) 



Differentiating eq. HA5|I on parameter /?' and using eqs. 0, ® we have 

^ (A6) 



— — = kZ 1 



On the other hand differentiating eq. iQJ on /?' and using distribution function ljA3|l and eqs. © we have 



dS kZ— " 



d(5' l-q 



-^^{l-{l-q)l3'E)+E 



(A7) 
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where was used dEi/djS' = 0. Comparing eqs. ljA6|l and l|A7|l the foUowing equahty is vahd 

ds dE 

Then the parameter (3' is expressed from the temperature 

1 _dS dS/d[3' , 



and 



:i _qZ 1 



kT 



(A8) 



(A9) 



(AlO) 



In this respect it should be noted that in the paper of Q.A. Wang multipHer q from eq. ljA10|) was lost as 
the entropy (|A5(I was directly differentiated from E without taking into account that the energy E is not an 
independent variable of canonical ensemble (T, V, N) . 

Substituting eq. (|A10|) into eq. IjASp and introducing the function A in the following form (see Section II) 



A EE kTq-^ "—^ =E- q-^TS 

1 - q 

for distribution function l)A3|l and norm equation (jA4|) we have 



Pi 



kTq- 



(All) 



(A12) 



and 



E 



1 + (1 - g) 



A-E,, 



kTq- 



(A13) 



Thus the minimization procedure of two different functionals H14|l and I|A1|I results in equivalency of distribution 
functions (O and HA12|) . 

APPENDIX B: THE CLASSICAL IDEAL GAS IN WANG'S REPRESENTATION 

1. Case g > 1 

In the case 5 > 1 it is convenient to use the following integral representation of the Euler's Gamma func- 
tion [li, 



QO 

[ ty-^e-^'^dt, Rex > 0, Rey > 0. 

r(y) J 



(Bl) 



To solve eqs. (|32|l one should substitute the values of a; = 1 + (1 — q)qPrj ^{E ~ Ei) and y = q/{q — 1) in 

the case of eq. ()31|) and the values of x from the above and y = 1/(^—1) in the case of eq. (I32() into formula 
IbH . We get 







(B2) 



where (3 — 1/kT and Zi,n is the canonical partition function of the ideal Boltzmann gas of N identical particles 
with the argument (3 t{q — l)q(3ri~^ 



Zi,N{t{q - l)qPr^) = {t{q - l)q?7"')"^'^^l,iv(/3). 



(B3) 
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Taking into account ljBl|) and ljB3p . one can carry out the integration of integrand ljB2|l with two values of y 
from the eqs. igU, 



- IN) 



[l + {l-q)q(}r^-^E\ —+'! 



N 



(g(g-l))i^r(^) 

Solving the eqs. (|B4|I . (|B5(I relative to the unknown variables E and 77, we have the following expressions 



E = -7V(g/3)-S, 



?7 = 



l + [l-q)-N 



- IN) 

(g(g-l))i^r(-^) 



(B4) 
(B5) 

(B6) 
(B7) 



which are vahd under the conditions -3N/2 + q/(q - 1) > and -3N/2 + l/{q-l) > 0. 

The one-particle distribution function are calculated from the eqs. (|78f) . (|80|) by substituting the values of 



x = 1 + (1 - q)qPr] ^{E - J2iLiPi f^^^) ^^'^ D — l/il ^ 1) "^to eq. 



f{p) = 



i^-l e-t[l+(l-,Wv-\E-^)] Z,^M-Mq-l)qf3r')dt. 



N{2^n)^ r(^) 

u 

Taking into account eqs. (|B1|I and (|B3|I . we get for the integrant (|B8() 



(B8) 



f{p) 



a \3/2 p,'_2 3/^_2)-) 



27rTO 



l + (l-g)g/377-i(i^-|^) 



(g(g-l))i(^-i)r(^) 

T^+i(^-i) 



(B9) 



Substituting the eqs. (jB6|) and ljB7p into expression ljB9|l we obtain the one-particle distribution function (|81|l 
with the coefficient ()82|l under restrictions described in the main text. 



2. Case q<l 



In the case q < 1 should be used the Hankel's contour integral in the complex plane [T3 making the trans- 
formation t ^ tx with X real and positive [TH . 



v-i 



2tt 



r(2/)T^ ^(-0"^e-*^di, i?ex > 0, |y| < 00. 



(BIO) 



To solve eqs. (j^ . (|^ one should substitute the values of x = 1 + (1 — q)ql3i]^^{E — Ei) and y = 1/(1 — 
in the case of eq. (PT|l and the values of x from the above and y — (2 — g)/(l — g) in the case of eq. into 
formula ljB10|l . We get 



=r(y)^ i(_i)-.e-*[i+(i-«)«/3')-^^^l Z,^M{t{q-l)qf3r^)dt, 



(Bll) 



After calculations similar to eqs. (|B2|l - ljB7|l for the energy E and the function -q, we have 



E - -N{qP)-^r^, 



l + {l-q)-N 



(g(i-g))-i^r(^) \-T^+f" 



(B12) 
(B13) 
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The one-particle distribution function are calculated from the eqs. (|78|) . (|80|) by substituting the values of 
X = 1 + (1 - q)q/3T]-^{E - E^iPf/2"i) and y = 1/(1 - q) into eq. |BT0)| 



fip) - N(Zhr '^^^^T^q^i fi-tr^ ^-tll+il-M-HE^^)] z,^^_,{t{q - l)qP^'')dt. (B14) 

c 



Taking into account eqs. (|B10|I and l|B3p . we get for the integrant l)B14(l 



fip) = 



V 2'K'm 



3/2 



l + {l-q)qPn~\E-^) 



(g(i-g))-i(^-i)r(^) 
r(T^ + 1(^-1)) 



(B15) 



Substituting the eqs. ljB12p and IjBlSp into expression (jB15|) we obtain the one-particle distribution function 
(ISU with the coefhcient 
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FIG. 1: The one-particle distribution function for classical ideal gas of A'^ = 500 identical nucleons in extensive Renyi 
statistics at the temperature T = 100 MeV for different values of q: q = 1 (the second line from below) correspond to 
the Maxwell's distribution function, other lines to the Renyi distribution function. 
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FIG. 2: The one-particle distribution functions for classical ideal gas of A'^ = 504, 500, 450, 400 and 1 nucleons from the 
top down in extensive Renyi statistics at the temperature T = 100 MeV and q = 0.99868. The value of TV is restricted 
to iV < I i_g~504 38 at this choice of the Tsallis index. 



